home *** CD-ROM | disk | FTP | other *** search
/ Just Call Me Internet / Just Call Me Internet.iso / com / computer / casio_st / prog_fx / divers / simplex.sha < prev    next >
Encoding:
Text File  |  1994-03-01  |  1.6 KB  |  72 lines
  1. 10 REM METHODE DU SIMPLEX
  2. 20 REM (C) WILLYTEL JUIN 89
  3. 30 REM 24H/7H AU 62.93.74.05
  4. 40 REM ENTREE VALEURS
  5. 50 A$="":INPUT "CLEAR ALL ?";A$:IF A$="O" THEN CLEAR
  6. 60 INPUT "Nombre de var:";NVAR
  7. 70 INPUT "Nombre d'equation:";NEQ
  8. 80 XX=NVAR+3
  9. 90 YY=NEQ+1
  10. 100 DIM T(XX,YY)
  11. 110 DIM R(YY)
  12. 120 REM ENTREE TABLEAU
  13. 130 FOR Q=1 TO YY-1
  14. 140 FOR P=3 TO XX-1
  15. 150 PRINT P-2;" ";Q;CHR$(13);:INPUT "Coef eq:";T(P,Q)
  16. 160 NEXT P
  17. 170 INPUT "Resultat:";T(P,Q)
  18. 180 NEXT Q
  19. 190 REM ENTREE Z
  20. 200 FOR P=3 TO XX-1
  21. 210 PRINT P-2;CHR$(13);:INPUT "Coef Z:";T(P,YY)
  22. 220 NEXT P
  23. 230 REM ENTREE SOLUTION BASE
  24. 240 FOR P=1 TO YY-1
  25. 250 INPUT "X(i)?:";T(2,P)
  26. 260 T(1,P)=T(T(2,P)+2,YY)
  27. 270 T(XX,YY)=T(XX,YY)+T(1,P)*T(XX,P)
  28. 280 NEXT P
  29. 290 GOSUB 630
  30. 300 REM CALCUL Z MAX
  31. 310 MZ=0
  32. 320 FOR P=3 TO XX-1
  33. 330 IF T(P,YY)>MZ THEN MZ=T(P,YY):J=P
  34. 340 NEXT P
  35. 350 PRINT "Vecteur entre:";J-2
  36. 360 REM CAL VECT SORT
  37. 370 MIN=9*10^30
  38. 380 FOR P=1 TO YY-1
  39. 390 IF T(J,P)=0 THEN R(P)=9*10^30
  40. 400 IF T(J,P)<>0 THEN R(P)=T(XX,P)/T(J,P)
  41. 410 IF R(P)>0 AND R(P)<MIN THEN MIN=R(P):L1=P
  42. 420 NEXT P
  43. 430 PRINT "Vecteur sortant:";T(2,L1)
  44. 440 T(2,L1)=J-2
  45. 450 T(1,L1)=T(T(2,L1)+2,YY)
  46. 460 PAS1=T(J,L1)
  47. 470 FOR P=3 TO XX
  48. 480 T(P,L1)=T(P,L1)/PAS1
  49. 490 NEXT P
  50. 500 FOR P=1 TO YY-1
  51. 510 IF T(J,P)=0 OR P=L1 THEN GOTO 520
  52. 520 CO=T(J,P)
  53. 530 FOR Q=3 TO XX
  54. 540 T(Q,P)=T(Q,P)-CO*T(Q,L1)
  55. 550 NEXT Q
  56. 560 NEXT P
  57. 570 CC=T(J,YY)
  58. 580 FOR P=3 TO XX-1
  59. 590 T(P,YY)=T(P,YY)-T(P,L1)*CC
  60. 600 NEXT P
  61. 610 T(XX,YY)=T(XX,YY)+T(XX,L1)*CC
  62. 620 GOSUB 630:GOTO 300
  63. 630 REM AFFICHAGE
  64. 640 FOR P=1 TO YY
  65. 650 FOR Q=1 TO XX
  66. 660 PRINT T(Q,P);" ";
  67. 670 NEXT Q
  68. 680 PRINT" "
  69. 690 NEXT P
  70. 700 INPUT "Encore?";A$
  71. 710 RETURN
  72.